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Abstract 

Fermion systems with more than two components can exhibit pairing condensates of much more 
complex structure than the well-known single BCS condensate of spin-up and spin-down fermions. 
In the framework of the exactly solvable S0{8) Richardson-Gaudin (RG) model with S'C/(4)- 
symmetric Hamiltonians, we show that the BCS approximation remains valid in the thermodynamic 
limit of large systems for describing the ground state energy and the canonical and quasiparticle 
excitation gaps. Correlations beyond BCS pairing give rise to a spectrum of collective excitations, 
but these do not affect the bulk energy and quasiparticle gaps. 
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I. INTRODUCTION 



3, la I wit 



nuclei 



The exact solutions of the pairing Hamiltonian with complete orbital degeneracy are 
available as simple analytic formulas not only for ordinary 5 = pairing of spin 1/2 particles 
with SU{2) symmetry [1], but also for systems of particles having additional degrees of 
reedom with respect to an internal quantum number, e.g. color degeneracy in quark models 
uh SO{Q) symmetry, or the isospin degeneracy of neutrons and protons in atomic 
sl with 50(8) symmetry. To fully describe pairing in nuclei, we should consider an 
isospin-invariant pairing Hamiltonian made of the three 5 = pair creation and annihilation 
operators transforming under isospin as T = 1 closing an 50(5) algebra. We call this the 
isovector pairing Hamiltonian. In nuclear systems where the number of neutrons is close to 
the number of protons Z, we should also consider the T = pairing leading to a Hamiltonian 
made of the six pair operators that carry 5 = 0,T = 1 and 5 = 1,T = 0. We show that 
the bulk limit of the pair correlation energy in the systems with neutrons and protons is 
equal to sum of the pairing energies of the two kinds of nucleons in isolation. Also, the gap 
in quasiparticle excitation energies is the same. This demonstration supports the use of the 
BCS theory to calculate pairing energies for the bulk limit of fermion systems with more 
than two degenerate degrees of freedom. 

The exactly solvable pairing Hamiltonians are all based on interactions that have a simple 
operator structure. The first we consider is the pure seniority 50(8) Hamiltonian of T = 0, 1 
pairing, for which analytic expressions for the energy can be derived. This model has 
been applied to nuclear physics in ref. g]. More realistic but still exactly solvable by 
algebraic methods are the RG Hamiltonians with include non-degenerate single particle 
levels (for a review see m). More recently, these exactly solvable models have been extended 
to higher rank algebras Is , 10 1 describing T = 1 pairing, color pairing and T = 0, 1 pairing 
respectively. For the case of T = 0, 1 we examine both exactly solvable sets, the degenerate 
orbital model and the RG model with non-degenerate single particle orbits, in turn. 
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II. THE MULTICOMPONENT SENIORITY HAMILTONIAN 



The pairing Hamiltonian with degenerate orbitals can be solved analytically if the inter- 
action has a separable structure with pair operators. The Hamiltonian is given by 

n 

H = -GY,Y.Pl^P-r (1) 

Here VL is the number of spatial orbitals and is a pair creation operator. The label a indi- 
cates the quantum numbers associated with the fermion internal degrees of freedom. More 
specifically, for the nuclear physics application, we distinguish isospin quantum numbers of 
the pairing operators and write the Hamiltonian as 

h = -gY, Y.^pIp,^ + Q\^,3). (2) 

with 

In the expression for the pair creation operators P\Q\ the first column in the coupling 
refers to spin and the second to isospin. The six pair creation and annihilation operators 
together with the spin operator, the isospin operator and the spin-isospin tensor close the 
5*0(8) algebra. The Hamiltonian ([2]) with equal strength for the two pairings (T = and 
T = 1) is invariant under S'f/(4) transformations of the fermion basis. 

Moreover, the proton-neutron pairing Hamiltonian ([2]) is completely equivalent to a spin 



= 3/2 fermion pairing 



hyperfine components 



iamiltonian appropriate to describe trapped cold atoms with four 



The pair operators that close the S0(8) algebra in this case are 
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where the spin 3/2 fermion pairs couple to total spin 5 = (monopole S pair) and 5 = 2 
(quadrupole D pair), and the SU{A) symmetric Hamiltonian reads 

n 2 

H = -GY,{S\S,+ Y,dId,3)- (4) 

ij n=-2 

There is a one-to-one correspondence between Hamiltonians (12]) and (jlj) if the pair op- 
erators are related as QI_i = -D±2) P±i = -^±1; -^o — ^'h Qo — ■ Therefore, all 
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conclusions we will extract from the proton-neutron S'f/(4) symmetric pairing Hamiltoni- 
ans are equally valid for the monopole and quadrupole pairing Hamiltonians of spin 3/2 
fermions. 

Coming back to the Hamiltonian ([2]), we shall also consider the identical-particle pairing 
Hamiltonian keeping only the two P-type operators P|^,p1^^, and the isovector pairing 
Hamiltonian keeping all three P-type operators. To define a bulk limit for Hamiltonian of 
the form Eq. (1), we let both the number of spatial orbitals Q and the number of particles 

go to infinity, but keeping the fractional filling fixed at some value /. The bulk limit also 
requires that the energy scale with the size of the system. This implies that G vary inversely 
to f2 for large values of f2, i.e. 

with g a constant. 



A. Identical-particle pairing 

As a warm-up exercise, we begin with the well-known case of the seniority Hamiltonian 
for the two-component fermion system, eg. only neutrons or only protons. The formula for 
the energy of the identical-particle seniority Hamiltonian is jl|, Eq. 6.12] 

E,{Nr) = -^G{Nr-v){2n-Nr-v + 2), (6) 

where v is the seniority and Nr = Np, Nn is the number of nucleons of isospin r. The ground 
state has seniority v = 0. In terms of g and / = Nt-/2Q, the energy per particle can be 
expressed 

^° ^fl-/ + AV (7) 



iV, 2^ V ^/ 
The corresponding expressions for the energy in the BCS theory arefll, Eq. 6.62] 

EBcs = -\GN(^2n-N+^^ (8) 

and 

^ = -^^7(1-/ + //^)- (9) 

The energies in eq. ([7]) and ([8]) only differ by a term of order 0{g/Q). Thus, the BCS theory 
for identical particles and the seniority Hamiltonian is exact in the bulk limit. 
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We now ask about the accuracy of Eq. ([8]) in the more general context where both 
neutrons and protons are together. For simphcity, we only consider the case of equal numbers 
of neutrons and protons, 

N 

N, = Nr, = -. (10) 

The naive extension of the BCS theory is to simply add the energies of the neutron and 
proton condensates independently. Then Eq. ([8]) would still be valid. 



B. Isovector and S'?7(4)-symmetric Hamiltonians 

The corresponding formulas for isovector pairing are given in ref. js], Eq. (9)]. The 
eigenvalues are 

E,^t{N) = -Ig{N -v){An-N-v + Q) + ^T{T + 1); (11) 
o Z 

where T is the isospin of the eigenf unction. Let us assume that A^^ and Np are even. Then 
the ground state has seniority v = and isospin T = 0. The energy per particle is 

^ = -l,fl-/+±V (12) 

One sees that the leading term is the same as in the identical-particle Hamiltonian. The 
finite-size correction is different, giving more binding for isovector pairing. 

The general formula for the energies of the S'f/(4)-symmetric Hamiltonian is derived in 
ref. [4]. It is rather complicated for general seniority, but for f = it reduces to a form very 
similar to Eq. ([7]) 

Eom{N) = -IcNim - iV + 12) + ^A2(A2 + 4). (13) 

o Z 

Here A2 is the spin-isospin label of of SU{A). In the ground state A2 = and the energy per 
particle is 

Again the bulk limit is the same as in identical-particle pairing. Note also that the finite-size 
correction is larger than either of the two other Hamiltonians. 
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C. Quasiparticle and collective excitations 



In the BCS theory of the identical-particle Hamiltonian, all excited states are quasi- 
particle excitations, whose energies can be calculated by blocking individual orbitals from 
participating in the condensate. The lowest excitation is a two-quasiparticle state; its energy 
is obtained by substituting fl ^ fl — 2 and N ^ N — 2 in Eq. ([8]). The resulting excitation 
energy 

E2g = g + 0{n'') (15) 

exhibits a finite gap in the excitation energy spectrum that persists to the thermodynamic 
limit. For the exact solutions, the two-quasiparticle excitations are given by the states of 
seniority v = 2. Formulas identical to Eq. f|T5l) are obtained by taking differences of energies 
between states of f = 2 and v = and keeping the other quantum numbers fixed. 

The other important manifestation of the BCS condensation is the number-parity depen- 
dence of the binding energies. The systems with odd are described as one-quasiparticle 
states in BCS and as states with seniority f = 1 in the exact solutions. In all cases the 
change in energy adding one particle is given by 

Eg = E{N + 1) - E{N) = I + 0{n-^) {N even) (16) 

= -| + oisr^) {N odd) 

The pairing gap A is conventionally defined as half the second-order difference of binding 
energies, with a sign chosen to have A a positive quantity. The formula is 

Ao-e{N) = ^ {2E{N) - E{N - 1) - E{N + 1)) (17) 

for odd. Again, the BCS and the exact solutions all have the same thermodynamic limit, 

A,_, = | + 0(r]-i). (18) 

While the quasiparticle excitations are the same in BCS and the exactly solvable Hamilto- 
nians, the latter have additional excitations corresponding to the other quantum numbers. 
For nonzero values of T in the isovector pairing Hamiltonian, the lowest excitation has 
V = 0,T = 2 with an excitation energy of 

£^0,2 — £^0,0 = (19) 
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This is inside the quasiparticle gap and goes to zero in the bulk hmit. The behavior of the 
5'f/(4)-invariant Hamiltonian is very similar; the excitations with v = and finite A2 are 
inside the quasiparticle gap and have vanishing excitation energy in the bulk limit. 

The collective excitations can be viewed as arising from the degeneracy of the BCS so- 
lutions with respect to choice of the paired orbitals. For the isovector Hamiltonian, the 
pair condensate wave function may be rotated in isospin space to give a degenerate ground 
state; the wave function transforms as the [11] representation of the SU{2) of isospin group. 
Similarly, for the S'[/(4)-symmetric Hamiltonian, the condensate transforms as the [11] rep- 
resentation of that group. This is all contained in the manifold of solutions of the Hartree- 
Fock-Bogoliubov equations. These equations give the natural generalization of BCS to 
permit arbitrary pair assignments in constructing the condensate. 



III. THE SO{8) RICHARDSON-GAUDIN MODEL 

To make a more realistic model for pairing in large systems, the Hamiltonian needs to 
include a single-particle term as well as the interaction. Richardson has given a construc- 
tion of the exact many-body wave function for two-component fermion systems with an 
arbitrary single-particle Hamiltonian and the separable pairing form for the interaction 
The combination of the Richardson solution with the integrable Gaudin magnet model 
gave rise to a large family of fully integrable and exactly solvable models based on the 
SU{2) algebra named Richardson-Gaudin models 7|. More recently the RG models have 
been generalized to simple algebras of arbitrary rank. In particular, the 5*0(8) RG models 
describe isoscalar-isovector pairing Hamiltonians with non-degenerate single particle orbits 
[lo| . In the simplest version of the model the Hamiltonian is S'f/(4)-symmetric, however 
symmetry breaking terms can be added to the single particle Hamiltonian within the exact 
solution. The exact solution of the S0{8) model is obtained by solving four sets of algebraic 
nonlinear coupled equations in terms of four sets of unknown spectral parameters. Each 
independent solution of the coupled set of equations completely determines an eigenstate 
and the corresponding eigenvalue. One of the sets of spectral parameters is composed of the 
pair energies that fix the structure of the correlated pairs in the single particle basis as in 
the SU{2) case. The other three sets of spectral parameters are responsible for the struc- 
ture of the eigenstates in the space of the internal degrees of freedom (spin-isospin). The 
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thermodynamic limit of the exact Richardson SU{2) solution has been obtained making use 
of an exact mapping between t he q uantum many-body problem and a classical electrostatic 



problem in two dimensions [15[ . This mapping is difficult to pursue for higher rank RG 
models due to the proliferation of different species of charges (spectral parameters). Instead, 
we will extrapolate to the bulk limit the numerical exact results that can be obtained for 
quite large systems. 

The SU{A) symmetric RG pairing Hamiltonian that we consider is 

n n 

H = Y.^^N^-|^ Y.{PlP,, + Dl^D,,), (20) 

i ijfi 

where iVj is the number operator counting all particles in level i and the isovector P^- and 
isoscalar Q^^^ are defined in ([3]). The single particle energies are taken equally spaced as 
Ei = ^{i — 1) with i an integer in the interval We carried out the calculation for 

quarter filling (A^ = Q) and interaction strength g = 0.15, considering systems of different 
sizes jl6!]. For these system parameters the fraction of the condensate is ~ 0.54. This value 
is estimated by counting the ratio of complex pair energies from the exact solution to the 



total number of pairs. As discussed in 



15| a complex pair energy implies the formation of 



a correlated Cooper pair, while the rest of pair energies, being real, describe free fermions. 
A. The BCS limit 

In the thermodynamic limit the BCS equations for the Hamiltonian Eq. fl20|) reduce to 

e-X 




+ {e- Xy 



I de 



1 1 

— n.F = —. 

v/A2 + {e- Xy 9 

and for the BCS energy we have : 




N 

The solution for the given parameter values has chemical potential A = 0.12468144 and 
jap parameter A = 0.015466976. The quasiparticle energy of level q in BCS is given by 



Eg = J{e,-Xy + A^ + X. (21) 
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TABLE I: Ground-state energy, quasi-particle energy, and different gaps as defined in the text for 
the 5'C/(4)-symmetric pairing Hamiltonian in the cubic expansion, with a comparison to bulk 
BCS Umit. 





Method 


a 


b 


c 


d 


E 
N 


Exact 
BCS 


0.062022149 
0.062022154 


-0.597581 


1.278831 


-11.1571 




Exact 
BCS 


0.140148 
0.140151 


-0.479740 


-10.0327 


-1107.25 




Exact 
BCS 


0.0154637 
0.0154669 


-0.699890 


-2.24642 


-1066.63 




Exact 
BCS 


0.0154672 
0.0154669 


0.0961964 


2.59458 


-257.910 



B. Thermodynamic limit of the exact solution 

As before, we want to compare the exact solution to the BCS approximation in the 
thermodynamic limit, N,Q ^ oo with N/Q constant. We first focus on the ground-state 
energy per-particle. We assume it is linear in in the thermodynamic limit, and we expand 

E{N)/N as 

E. = a + b/N + c/N^ + d/N^ + O (l/N^) . (22) 

For our analysis we find the exact numerical solutions for A^ in the range 160 < N < 1000 
and fit the parameters of the expansion Eq. (122|) . For details of how we solve the equations 
for the 50(8) RG model, see ref. [l^. The results are shown in Table I and Fig. 1. We 
see that the fit with a cubic polynomial reproduces the BCS energy with seven significant 
figures. We next examine the quasiparticle energy Eg, calculated as the binding difference 
between the systems with An and 4n + 1 particles, 

Eg{An) = E{An + 1) - E{An) ^a + h/N + c/N'^ + d/N^, (23) 

We show in the Table the 1/A^ expansion for this quantity. In our case we choose the 
blocked orbital as the lowest unoccupied level corresponding to = 1/8. We expect that 
the constant term will be the same for addition and removal energies, but the higher order 
terms in the expansion will differ. 
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FIG. 1: Ground-sate energy per-particle of the SC/(4)-symmetric fermion systems as a function of 
particle number A^. Black squares show the calculated results and the solid line is the cubic fit of 
Eq. (j22l) . The open circle is the bulk BCS limit. 

We also report in the Table the pairing gap computed in two different ways. The first is 
the binding energy difference Ao_e center at particle number = 4n + 1, 



where tt-j^t is the occupation probability of a fermion with spin a and isospin r in level 
i. Note that this definition coincides with the BCS gap if the occupation probabilities are 
calculated in the BCS approximation. 

The cubic fits to the quasiparticle energy and the odd-even and canonical gaps are shown 
in Table 1. The precision with which the thermodynamic limit of the exact solution ap- 
proaches the BCS results for these quantities is slightly lower than that of the ground state 
energy. The reason is that the quasiparticle energy and the odd-even gap result from the 
difference of big numbers (ground state energies of finite systems) while the canonical gap 
involves a further calculation of the occupation numbers once the solution of the spectral 




(24) 



Another, more theoretical definition is the canonical gap, 




(25) 
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FIG. 2: Gaps of the 5C/(4)-symmetric fermion systems as a function of particle number N. The 
continuous hues are the cubic fits. The empty diamonds is the bulk BCS gap parameter. 

parameters is obtained. In any case the precision obtained, five or six significant figures, is 
enough to guarantee that the thermodynamic limit of the exact RG model coincides with 
the bulk BCS values for the four magnitudes studied. The convergence of the two gaps 
with the size of the system is shown in Fig. 2. One can see that for these quantities the 
expansion is rapidly convergent as for the ground state energy, and that it can be reliably 
extrapolated. A noteworthy fact is that both definitions of the gap, the odd-even energy 
difference (l24l) and the canonical gap ( |25l) . converge to the bulk BCS gap, in spite of the 
fact that the approaching direction of convergence is opposite. 

IV. PERSPECTIVE 

The excellent agreement between the BCS and the extrapolated exact solutions of our 
numerical example shows that higher order correlations on top of a single pair condensate 
are not important for the bulk ground state characteristics of many-component fermion 
systems, when the Hamiltonian is of the pairing type given in Eq. (l20l) . However, the 
interaction in that Hamiltonian is still rather restrictive. For example, one could imagine a 
two-particle interaction of the 4-component system that supported a 4-particle bound state 
(the alpha particle in nuclear physics) but no heavier ones. Then the many-body system in 
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the Bose-Einstein condensate limit would be a gas of the 4-particle clusters. Analogously, in 
the weak coupling regime there might be stronger quartet correlations than in the current 
RG models. It might be interesting to ask what conditions would be needed to give the 
two-particle interaction a capability to produce extensive higher-order correlation effects. 

Although for pairing Hamiltonians like (1201) the higher correlations do not affect the bulk 
gaps and per-particle condensation energy, there are still collective excitations present in the 
system than are beyond the reach of the BCS approximation. This is seen most clearly in 
the analytic formulas of the pure seniority Hamiltonians discussed in Sect. II. It might be in- 
teresting to see how well the standard methods of many-body theory perform for calculating 
the effects of the collective degrees of freedom starting from the HFB solutions. For example, 
collective effects are often treated in nuclear theory by the Quasi-Particle Random Phase 
Approximation or by the Generator Coordinate Method. It fact, it has been shown already 
that for the 2-component fermion system the particle-particle RPA and QRPA give very 
accurate corrections to BCS for the total ground-state energy [17|, ll8|] of finite systems in 
the fluctuation dominated regime and in the superconducting regime respectively. However, 
close to the phase transition region both approaches cannot reproduce the large quantum 
fluctuations. It would be interesting to carry out these analysis in the thermodynamic limit, 
as well as to consider quantum corrections beyond BCS to multicomponent fermion systems 
having the exact solution of the higher rank RG models as a benchmark reference. 
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